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The first order phase transition ordinary occurs in the systems with many 
different condensating components. During the process of nucleation, i.e. 
formation of droplets, the heat of condensation is extracted which changes 
the rate of nucleation. The theoretical description has to take into account 
these two features. So, one can see that creation of the nonisothermal theory 
for the multicomponent nucleation is rather actual. 

The adequate theoretical description of the nucleation stage was given for 
many different situations starting from the pure isothermal nucleation of one 
component (one substance). Even in this situation described by the so-called 
classical theory of nucleation there is no coincidence between theoretical pre- 
dictions and experimental results. Nevertheless, this disagreement can not 
be the reason to reject all further modifications of the theory to grasp the ef- 
fects of the heat extraction and to extend the theory for the multicomponent 
case. 

The first theory where the thermal effects were taken into account rather 
adequately was presented by Kantrowitz In this publication embryos 
of given size were characterized by the unique average temperature. Only 
the so-called "weak thermal effects" were considered there. But the really 
essential decrease of the nucleation rate occurs under the "strong thermal 
effects" considered by Kuni 0. 

Certainly, there exists the distribution of the embryos of given size over 
the temperature. Feder et al. took this phenomena into account for the 
weak thermal effects. The energy distribution of embryos under the strong 
thermal effects was investigated by Kuni and Grinin ||^. 
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The standard way to construct kinetic equation is to use the Fokker- 
Planck approximation. Then kinetic equation can be solved in the stationary 
situation by approach of Langer [^] or by approach of Kuni et al. [^] in the 
nonstationary situation. 

The methods described in allow to study some different situations. 
Namely the situation of the binary situation was studded by Lazaridis and 
Drossinos under the Fokker-Planck approximation. 

Concerning the mentioned publication one has to note that due to a 
rather big quantity of the molecular heat extracted in the act of condensation 
the restriction by the Fokker-Planck approximation isn't sufficient. Then 
kinetic equation contains high derivatives and one has to use the Chapman- 
Enskog procedure to solve it. 

Until the current moment one can say that the most advanced approach 
for the nonisothermal nucleation is the application of the Chapman-Enskog 
procedure which was suggested by Kuni, Grinin in the case of one component 
nucleation 0]. In |P] this approach was spread to the binary nucleation. 
Unfortunately it was not done in the appropriate way. So, the subject of 
this publication appears rather naturally in the context of the theoretical 
methods development. 

Particularly, this publication is aimed to present the nonisothermal theory 
of nucleation in the multicomponent mixture and to correct the errors in 
We shall present the self consistent theory which gives the analytical 
expression for the nucleation rate. Some parts of the theory where it is quite 
analogous to the known results (for example, the solution of kinetic equation 
in the Fokker-Planck approximation) are omitted. 

In this publication we shall use all standard definitions of the classical 
nucleation theory, consider the unit volume, take all values with the en- 
ergy dimension in the units of the elementary thermal energy kbT [kh is the 
Boltzmann constant and T is the temperature), and all values with the heat 
capacity dimension in the kf, units. 

1 Substance exchange 

Here we shall derive kinetic equation from the balance equation for the dis- 
tribution n over the embryos sizes, this balance equation can be written in 
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the following form 

where V is the operator associated with the substance exchange between 
the embryo and the environment (one also has to take into account here the 
energy exchange due to the extraction of the condensation heat) and G is 
the operator associated with the energy exchange due to the difference of the 
temperature T of the embryo from the temperature Tq of the environment 
(the effects of the heat extraction are already included into V). 

The function n is the number of embryos of given size in a unit volume. 
Later we shall specify variables of the embryos state description. 

At first we shall study kinetic equation without thermal relaxation by a 
passive gas, i.e. in the form 

and then we shall study thermal relaxation, i.e. operator G. 

Suppose that there exists a vapor mixture with condensating com- 
ponents. All these components can be found both in the vapor and in the 
embryo (in contrary of the passive gases which can be found only in the vapor 
phase). The distribution n{{i'i},t) of embryos is the function of io variables 
of the number of molecules z/j of the ith component (substance) inside the 
embryo and also the function of time t. 

The process of the absorption of a molecule of jth component occurs with 
intensity W^aj cond where Wj' is the intensity of collision of the given embryo 
with an arbitrary molecule of jth component and aj cond is the condensation 
coefficient for jth component. The most natural is to suppose that the length 
of the free motion in the gas media strongly exceeds the linear size of the 
embryo (so-called free molecular regime of the substance exchange). Then 
the value can be easily found by the gas kinetic theory 

W+ ~ Srij 

where S is the surface square of the embryo and nj is the molecular num- 
ber density of the jth component^. Due to the small relative size of the 

^One has also to note that in the nearcritical region all embryos are in the quasiequi- 
librium with the surrounding vapor. This quasiequilibrium also leads to the small unho- 
mogenities in the vapor phase. These unhomogeneties would lead to the diffusion flows on 
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nearcritical region^ we can see that Wj' remains practically constant in the 
nearcritical region. 

Any reliable information about the condensation coefficient aj cond is ab- 
sent (in the literature one can found rather different estimates). But we 
believe that aj cond is rather smooth function of the embryo state and can be 
regarded in the nearcritical region as some constant value^. 

The process of absorption of the molecule leads to the variation of the 
number of the molecules inside the embryo 

{M + 1} 

and also to the extraction of the condensation heat (3j (which is measured in 
the natural thermal units) which is going to increase the temperature of the 
embryo T by the value 

where Cj are the heat capacities per one molecule in the liquid phase (taken 
in units of fcf,), the sum is taken over all components. 

The value Wj~ depends only on the state of a vapor-gas mixture. Contrary 
to Wj' the intensity of the ejection of the molecule of the j'th component 
strongly depends on the temperature of the embryo. The embryo has to be 
characterized by the temperature T of the embryo^ or some function of the 
temperature. Instead of T one can introduce the value of additional energy 
E according to 

Jo 

where Tq is the temperature of the media, q are molecular specific heat^ 
expressed in units of kh- The evident advantage of E is that the equilibrium 
value coincides with the zero point. Now we shall normalize E in order to 

the embryo. That's why one can approximately say that the vapor around the nearcritical 
embryo is unperturbed and calculate the flow by the previous formula. 

^One can use the standard classical nucleation theory for this estimate. 

•^This follows also from the model where aj cond appears as the probability to overcome 
the energy barrier near the surface of the liquid phase. 

"^We suppose that the temperature relaxation inside occurs very rapidly. It can be 
justified by estimates analogous to 

^They are close to those defined under the constant pressure. 
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have no coefficient in the square form of the equihbrium distribution^ 
along additive energy near E = 0: 

rf ~ exp(— /i^) 

To find yU we shall start with Clapeyron-Clausius formula for the molecular 
number density rioo j of the saturated vapor under a plane surface of liquid 

T — T 

'^oo jiT) = J (To) exp(/3j ° ) 

J-o 

As far asQ for the planar surface 

Wf{noo j) = Wj'{noo i) ~ noo j 

one can come to| 

W-{T) _ ,{T) _ .^ T-Tq 
Wr{To) n^,{To) """P^^^' To ' 

One the other hand 

Wr{E)n%E) = W^n\E - pj)rf{E - 

and 

W-{E = 0)n%E = 0) = W+n'{-pj)n'{-pj) 

which leads to 

-^(-' ^ 

where F is the free energy of the embryos formation and it is taken into 
account that the equilibrium distribution n'^ ~ exp(— F). 

^We shall omit the normalizing factor tt^^/^ of the equilibrium distribution. It can be 
easily reconstructed anywhere. 

^/here we suppose that the surface tension doesn't depend on the embryos temperature. 
Details of this approximation can be found in |^ . 

^More carefully one has to go at first from the critical embryo to the plane surface, then 
use the mentioned relation for the plane surface and finally return to the critical embryo 
back. 
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Both approaches will coincide when 
ThenQ 



E 

/i 



One can write kinetic equation (|I|) in variables in the following 

form 



Ol = 1^ W^niiUi^j, Uj - 1}, - Tj) - 2^ Wj n{{vi^j, uj}, ^) 

j j 

j j 

where 



' iEiCiiyiV/' 
One can present the following split of the last equation 

— Q-^ — = i^yJAi^i^v ^3 - 1}. - ^3) - JAi^i^j^ ^3}^ N\ 

where the flow Jj is defined by 

JA{^^}^ /^) = W^n{{ui^j, Uj}, /i) - Wfn{{ui^j, + 1}, /i + tj) 

Now we have to substitute the finite differences by derivatives. 

One has to mention that the elementary steps 1 along z/j are small in com- 
parison with the characteristic scale corresponding to the essential variation 
of exponent of the free energy. This allows to substitute the finite difference 
along Ui only by the first derivative^ 

^More carefully it can be done in terms of the finite differences. 

^"^This produces certain restrictions which will limit later the Chapman-Enskog 
expansion 
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An elementary step tj along /i corresponds to essential violation of expo- 
nent of the free energy. So, one has to substitute the finite difference by the 
whole Tailor seria. As the result we have 

d 

We have to note that the possibility to substitute Jj({z/j^j, z/^ — l},/i) — 
Jjdiyi^j, I'j}, fi) by —■^Jji{i^i^j,i^j},fJ') can be made when we consider the 
situation near the quasistationary one. Then in one dimensional projection 
on Uj we shall get the small value for Jj({i/j-^j, z/^ — l},/i) — Jj ({i/j^j, 
which allows to substitute it only by the first derivative. 

The flow Jj can be expressed with the help of a function 

as 

The analogous substitution of the finite differences by the Taylor seria gives 

Jj{{iyi^j,iyj},fi) = iy+n'^({z/i},/i,t) 

To be close to the standard form of one component nonisothermal theory 
[^] we shall use instead of function / the following function 



exp(— /i^) 

The free energy F of the embryos formation can be split as 
F({z/},/i)=F({z/},/i = 0)+/i2 
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which gives 

fi{u}, /i) = P{{u}, fi) exp(F({z/}, ^ = 0)) 



One can present expression for Jj in terms of function P as 

d_ 



or 



exp(F({z/i i^j, + 1}, /X = 0) ^ 'ir\'Fh^^^^^^^i' + 

m=l '^A'" 



d 

Jji{'yi^j,'yj},f^) = W+n^{{ui}, fi)[-— exp{F{{u}, = 0)P{{ui^j,Uj}, fi) - 



Having introduced an operator 

one can present the last expression for Jj as 

Jj{{ui^j,Uj},fj,) = exp{-n^)LjP{{iyi^j,Uj},fi) - 

According to the smooth dependence along Vj one can substitute the 
argument Vj ± 1 of function|^ P by Vy 
Then as far as 

dF Id'^F l,dF^ _ dF Id'^F I ,dF ^ 



Later it will be seen that d/dvj produces some small parameter as far as the action of 
dF/dvj. We write this equation in the first two orders of this parameter. This corresponds 
to the order essential in the isothermal version of the theory. 
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then kinetic equation can be presented as 

exp(-/x^) 



EE 

3 ;=i 



dt 



oo ^ y Qi 



dF Id'^F 1,5F,2, , 2^..r+ ^ d"" r.^r 

E^[exp(V)L,P({..)..) - (1 + + ;^)-P(-M^)H^- E i;^,a^nM)\ 

The action of on exp(— /z^)-!/) where -0 is the arbitrary function is 

obviously given by 

Ql gl 

^exp(-/x2)V^ = exp{-id^){—^ -2/x)V 

Certainly one can not take 2/x away from ^ and has to consider — 
as sequential action of operators in brackets. This turn kinetic equation to 

dt ^! ^^/^ ' 

The last equation present the final form of kinetic equation. 

Wc decompose the finite differences along Vj until the second derivatives 
(or the second order of the small parameter) because in the classical theory 
of isothermal one-component nucleation two derivatives have to be taken into 
account (An account of the first derivative couldn't lead to the suitable rate 
of nucleation). 
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One can easily note that operators 



S -1 ^ 



OF Id'^F IdF. 



are absent in Really, in the first two steps of the Chapman- Enskog 
procedure described later these terms will not be essential. But it occurs only 
in frames of the Chapman-Enskog procedure and can not be seen directly 
from this equation. We shall call Sij, S2j, S^j as the shift operators. 

Now we shall turn to the thermal relaxation by the passive gas in order 
to include it in the presented equation. 



2 Thermal relaxation 

The physical reason to consider the interaction of the embryo with the passive 
gas is rather simple. Really, due to heat extraction the temperature of the 
embryo is higher than the temperature of environment and the embryo heats 
the molecules of the passive gas. Certainly, the temperature of the embryo 
falls which reduces an ejection rate. This has to be taken into account and 
the consideration of the interaction with the passive gas is important. 

In the previous consideration the condensation of the molecule can be 
described by some fixed values of the the condensation heat jSj. An obvi- 
ous restriction only by the regular term in the presence of the big quantities 
of the passive gas will lead to the thin spectrum in the energy scale of the 
(5— function form. Certainly, this doesn't coincide with the equilibrium distri- 
bution. Thus, one has to use at least the Fokker-Planck approximation. The 
physical reason is rather obvious - molecules of a passive have the different 
velocities and equilibrium distribution in energies. This has to be taken into 
account and leads at least to the Fokker-Planck approximation. As far as 
the variation of the energy in the elementary act of interaction is small in 
comparison with the characteristic scale of the variation of the equilibrium 
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distribution one can restrict this description by the Fokker-Planck approxi- 
mation. 

In Fokker-Planck approximation the kinetic equation can be written as 

dn d ^ d 
— — B—n—f 
dt Ofj, Ofi 

where B is kinetic coefficient. It can be determined by consideration of the 
limit situation where the last equation transforms into an equation only with 
the regular term 

dn ^ d (9 , ^ d ^ 

This form form has to be reproduced by the standard analysis. We begin 
with the balance equation 

where is the rate of coUisions of the given embryo with the molecules of 
the passive gas, the regular variation 5// is given by 

5/^ = 

and Cg is the molecule heat capacity of the passive gas. This leads to 

B^W+- 



2 E CjVj 

Now the kinetic equation can be written in the following form 

dn , Cg d ^ d 

— — W n — / 

dt 2Y^Cji'jdfi dn 

One has to put the coefficient of thermal accommodation aacc to in 
order to take into account that thermal accommodation occurs with some 
probability. 

The generalization of the previous equation on the mixture of passive 
gases leads to 
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where indexes with prime denote different passive gases. 

One has to talce into account that the condensating substances also take 
part in the coohng. With the probabihty (1 — CKc)cKacc the act of coohng takes 
place. The molecules accumulated by embryos also have to be taken into 
account. As the result one can get 

dn ^ ^ Cgi, d ^ d 

at ^ CjUj dji oji 

It is quite obvious that now to get the general kinetic equation we have to 

add the part associated with the condensating substance. So, it is necessary 
to formulate the part under consideration in terms of function P. Here one 
has to fulfill the same actions and get 

As the final result for the general kinetic equation one can get 

dP{{u,},f,) _ 
dt 

°° (—T Y r) °° r"* r)^ 
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Here indexes and j mark the different components of the condensating 
substances and index j' marks the different components of the passive sub- 
stances. 

One can easily note that the number of components of the condensating 
mixture doesn't act on the properties of passive gases. This hes in con- 
tradiction with results presented in [|| where the action of the passive gas is 
referred to the action of every component of condensating substance and then 
the direct summation over the condensation components is carried out. So, 
according to [|] one can speak about the separate cooling of different com- 
ponents (in kinetic sense, the droplet is cooling as a whole object). Here we 
speak about the common cooling of different components (in kinetic sense). 
The physical essence is another here. 



3 Estimates of operators 

Now we shall present the method to solve the last equation. The general 
scheme is well known. At first the extraction of the main operator with 
the well known eigenfunctions has to be presented. This main operator has 
to ensure the relaxation to the stationary state which allows to consider of 
the relaxation period. Such a structure allows later to apply the Chapman- 
Enskog procedure. 

An attempt to investigate the situation of the binary nonisothermal nu- 
cleation was made in [Q] but the initial kinetic equation was wrong. Contrary 
to [§ we shall use the correct kinetic equation. Here this equation is already 
generalized for the multicomponent case. 

Besides the new object of investigation the approach presented here has 
also some new principal features. 

One has to note the specific feature of the relaxation stage description 
fulfilled in the mentioned papers. When the main operator has only formal 
priority then a standard consideration of the relaxation stage requires the 
small value of specific parameter. This leads to the serious restriction of 
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the approach used in the mentioned papers. In the situations of intensive 
droplets formation this parameter isn't too small and the relaxation doesn't 
take place. Then the initial condition for the Chapman- Enskog procedure is 
violated. This doesn't allow to apply this procedure. 

We shall use another split on the r.h.s. of kinetic equation into the main 
operator and the additive one. As the result we come to some more compli- 
cate procedure with two sets of the main and additive operators. But still in 
such situation it will be possible to generalize the Chapman-Enskog proce- 
dure and come to the final formulas. It will be possible to get the relaxation 
to the stationary state without restriction used in [Q. One has to mention 
that the cited papers couldn't overcome the main nontrivial feature of the 
nonisothermal condensation - the main operator extracted in these papers 
has only formal priority based on the presence of factorials in denominators 
in the Taylor's expansion terms. So, one has to fulfill at least many steps in 
the Chapman-Enskog procedure. Here we shall present the way how to take 
into account the tails of these series and to come to the compact final results. 

Now it is worth mentioning the inclusion of the present analysis into the 
general scheme. 

One can note some specific features of the thermal effects in comparison 
with the general situation of the non Fokker-Planck evolution considered in 
[|T^. Namely these features allow us to go further in comparison with []10 



and to get the compact final formulas. 
These features are the following ones: 

• The temperature of the embryo can have an arbitrary value. 

• Non-Fokker-Planck evolution occurs along the temperature of the em- 
bryo. 



Non-Fokker-Planck evolution occurs under the constant value of /? 



r 



The third feature is rather important. Really, as far as we have the 
Clapeyron-Clausius relation we can reconstruct W~{T). Then on the base 
of (it is given by the simple gas kinetics formula) and W~ one can get 
the equilibrium distribution. The knowledge of the equilibrium distribution 
gives on the base of the Boltzmann formula the form of the free energy F 
of the embryo formation (the constant shift appeared from the normalizing 
factor of the equilibrium distribution isn't important). These constructions 
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result in a rather simple form of the free energy. In the arbitrary situation 
(see [jlO[ ) the form of the free energy can be more complicated and this causes 
the additional difficulties. 

The mentioned simple form of the free energy corresponds to the simplic- 
ity of transition from the function n to the function P defined by 

n = exp(— /i^)P 

Certainly exp(— /i^) represents here the equilibrium distribution and ap- 
pears due to the square character of the free energy. This leads to the P 
relaxation to a constant. To conserve such relaxation in the general more 
complicated situation one has to choose in the argument of exponent instead 
of /x^ another more complicated function which refiects the more complicated 
behavior of the free energy. As the result the eigenfunctions of the "main" 
operator will be unknown. 

In the general situation instead of fi in the combinations — 2/i + d/dn 
in the kinetic equation appear high powers of /x. This blocks the presented 
approach to get solution. 

We are going to act in frames of the macroscopic description of the free 
embryo. This leads to the big parameters 

Z/j c > 1 

for all components which are marked by index i. Index c corresponds to the 
critical embryo. 

The last inequality allows as it is shown in to state that 

• The Fokker-Planck approximation is valid to describe the evolution 
along z/j 

• The square approximation for the free energy along z/j in the near crit- 
ical region is valid. 

In [|T0[ all specific situations appear only when the derivative of the free 
energy along the concentration of the solution inside the embryo provides 
another big parameter. Certainly the essential values of this derivative are 
rather ordinary in the nature but one can not consider them as the big 
parameter going to infinity. 
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Really the derivative along Vi has a big value in comparison with the 



derivative along the steepens descent line [|T0[- This appears as the base 
for the hierarchy in the nearcritical region But the value of this big 

parameter isn't sufficient to compensate the influence of the big parameter 
z/j (as far as all Ui has one and the same power we shall drop the index i in 
the estimates). Namely, the halfwidth Az/ along u has the order i^^^, the 
halfwidth along the steepens descent line has the order i^l^^. But as far as 
Az/ is greater than 1 (it isn't so great as in the one component theory but it 
is still great) we can see that the differential form of the kinetic equation is 
valid. 

The last result can be directly seen from the explicit expression for the 
free energy of the embryos formation as the function of {z/j}. In the capillary 
approximation this expression can be written as 

j j 

Here Vj are the molecule volumes in the liquid phase, bj are the excesses of 
the chemical potentials, a is the renormalized surface tension. The surface 
of tension is put as to contain precisely the volume of the embryo. All Vj are 
supposed to have one and the same order, all Uc are also supposed to have 
one and the same order. 

One can easy note that the halfwidth Au ~ z/^/^ in the multicomponent 
theory differs from the same value in the one component theory Au ~ z/^/^. 
The reason is the interaction between components. This phenomena doesn't 
lie in contradiction with the general theory because the steepens descent line 
doesn't coincide with any Ui and the halfwidth along the steepens descent 
line coincides with the halfwidth along z/ in the one component theory. But 
it shows that the direct differentiation of F along z/j without the influence 
of the other components taken into account cannot lead to the really small 
parameter. An account of the mentioned interaction is rather difficult and 
it is more convenient to go to variables = z/j/ Y,j and k = c?!"^ ^j'^j- 
Then the direction along n coincides with the steepens descent line and due 
to the Gibbs-Duhem equation an account of the interaction is attained au- 
tomatically. 

In the set {^j}, k the form of the free energy is given by 

F ^ Bk- k^I^ 
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where B is some function of {^i}- The characteristic scale of n can be put as 
to coincide with the scale of z/j. 

To justify the validity of the square approximation one has to get the 
second and the third derivatives of the free energy. It is more convenient to 
use the last form of F. Then 



As far as the function B and it's derivatives don't contain any big parameter 
one can easy see that the action of the third term is small. That's why the 
square approximation for the behavior of F along is valid. The behavior of 
F along K is similar to the one-component case. The square approximation 
along K is, thus, valid. As the result the square approximation for F in the 
nearcritical region can be used. 

These results explain why the Fokker-Planck approximation is adopted 
for the description of the evolution along z/j. As for the evolution along 
temperature the Fokker-Planck approximation isn't sufficient. The reason is 
the existence of another big parameter. This parameter is 

Let's explain whyQ j3 can be regarded as the big parameter of the theory. 
When the temperature decreases from the value of the second order phase 
transition the value of of the molecule volume in the vapor phase grows 
and the value vi falls. So, far from the point of the second order phase 
transition^ one can come to 



^^The index will be omitted. 

^•^Namely in this situation some actual assumptions of the classical theory of nucleation 
are valid, for example the uncompressibility of liquid phase. 



d^F cPB 
and the halfwidth along is given by 



~ ^-'^\B"i{^,}))-'/' 




the action of the third term in the Taylor seria is given by 




Vv > Vi 
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The heat extracted in the phase transition can be presented as the dif- 
ference of entropies in two phases multiphed by temperature. Then one can 
use the standard representation of the entropy as the logarithm of the states 
number. The number of states^^ can be very approximately estimated as the 
volume occupied by the systemf^. Then the last strong inequality leads to 
the big value of /3. 

In reality one can not go very far from the second order phase transi- 
tion temperature because the new phase transition (crystallization) occurs. 
Nevertheless one has to say that (3 is the big parameter of the theory. 

Now we come to the direct solution of the kinetic equation. We can 
rearrange it in the following form 

^^i^ = D, + D. + D3 + D. (3) 
D, = -J2^L,Pi{u,},f^) 

] 1=1 " m=l (m^l) " 



EE' 

j 1=1 



d 



l\ 



2/i)%-5i,P({z/,},/i) 



a oo _m am 



Z3, = (4) 

^"^One can use quasiclassical approach. 

^^We can use very approximately the model of ideal gas. 
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i:|^(|-2.)'iy/s4|;,p({..},,) 

Now we shall estimate the actions of operators Di - D^. One can easy 
note that the differentiation along Ui can be estimated as 



dvi Avi 

which produces the small parameter 1/Az/j. One has to note that as has 
been already noted the value Ai/j differs from the analogous value in one 
component case which can lead to the error made in The values Ai/j don't 
estimate the size of the nearcritical region (The standard definition of the 
nearcritical region is given by inequality \F — Fc\ < 1. The infinite tails can 
be cut off to reduce the form of the nearcritical region to a rectangular one.). 
To estimate the size of the nearcritical region one can take derivatives of the 
free energy along k, ^j. As the result we come to the following convention: we 
use the notation l/Aj/j but keep in mind that the real small parameter will 
be I/Ak. Moreover one can not estimate the size of the nearcritical region 
by Az/j as it was done in but has to use the halfwidths along k and ^j. 

The value of derivative dF/dvi also contains the small parameter 1/Ai/j. 
So the action of Lj can be estimated as 

UP r^w; ^ 



AVi 



The differentiation along /i doesn't produce any small parameter. The 
characteristic value of /i is 1. 

Now we can calculate the powers of operators - D^. As the result we 
see that 

• Operator D4 is the smallest one. It has the order l/(Ai/j)^. 

• Operator D3 is small. It has the order 1/Ai/j. 

• Operators Di and D2 have one and the same order 1. Here operator 
Di has the formal priority because D2 has no terms without factorials 
in the denominators. 
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The main problem of the further analysis is that the main operator Di + 
D2 ensures relaxation to the state which gives zero flow of the embryos from 
the precritical to the postcritical region. So, this state leads to the zero 
value of the nucleation rate. To overcome this difficulty one can use the 
Chapman- Enskog procedure. 
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